Combination of random-barrier and random-trap models 

K. Mussawisade, T. Wichmann, and K.W. Kehr 
Institut fiir Festkorperforschung, Forschungszentrum Jiilich GmbH, 

52425 Jiilich, Germany 
(February 1, 2008) 

Abstract 

The temperature dependence of the diffusion coefficient of particles is stud- 
ied on lattices with disorder. A model is investigated with both trap and 
barrier disorder that was introduced before by Limoge and Bocquet (1990 
Phys. Rev. Lett. 65 60) to explain an Arrhenian temperature-dependence of 
the diffusion coefficient in amorphous substances. We have used a generalized 
effective-medium approximation (EMA) by introducing weighted transition 
rates as inferred from an exact expression for the diffusion coefficient in one- 
dimensional disordered chains. Monte Carlo simulations were made to check 
the validity of the approximations. Approximate Arrhenian behavior can 
be achieved in finite temperature intervals in three- and higher-dimensional 
lattices by adjusting the relative strengths of the barrier and trap disorder. 
Exact Arrhenian behavior of the diffusion coefficient can only be obtained in 

infinite dimensions. 
05.40.+j;05.60.+w;61.43.Fs 
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I. INTRODUCTION 



Many amorphous substances exhibit hnear behavior in an Arrhenius plot of the diffusion 
coefficient or the mobihty B of particles, where ln(D), or ln(i?), is presented as a function 
of the inverse temperature (3 = l/TkB- The observation of an Arrhenian temperature- 
dependence is not easily understandable from the theory of diffusion in disordered crystals. 
The commonly used models of diffusion of particles in lattices with disordered transition rates 
predict different behavior: The random-trap model predicts generally convex (downward) 
curvature of\n{D) versus (3, independent of the lattice dimension d, while the random-barrier 
model gives concave (upward) curvature for three- and higher-dimensional lattices. Limoge 
and Bocquet Jll suggested that the apparent Arrhenius behavior of \n{D) versus (3 might be 
due to a compensation of the effects of random barriers and random traps. The aim of this 
paper is the examination of this appealing proposition by analytical and numerical methods. 

The analytical arguments given by Limoge and Bocquet [|l| in support of a compensation 
of the effects of random barriers and of random traps are only partially satisfactory. They 
employed continuous-time random-walk theory and used a decoupling approximation when 
performing the disorder average. This procedure gives correct results for the random-trap 
model. For the random-barrier model, however, the ensuing results are not correct, as 
a consequence of the neglect of important backward correlations in the transitions of the 
particles. The failure of the decoupling approximation for the random-barrier model in 
comparison to the Monte Carlo results was already discused in [|l|. 

Recently a new exact result for the diffusion coefficient of particles in one- dimensional 
disordered chains has been found The expression contains transition rates that are 

weighted by the equilibrium occupancies of the sites; this is relevant for random site en- 
ergies, i. e., random traps. The insight obtained from this result can serve as a basis for 
approximate treatments of models with site and barrier disorder in higher dimensions, as 
will be shown here in Form of an effective-medium approximation (EMA). Very recently 
Limoge and Bocquet [^] tried to take the inherent backward correlations of the random- 
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barrier model into account. The differences between their and our results will be discussed 
below. 

In the remainder of the Introduction we present qualitative arguments for the different 
curvatures in the simple models. The downward curvature of ln{D) versus (3 for the random- 
trap model is easily understood: The deepest trap sites with the lowest energies dominate the 
behavior at the lowest temperatures. The convex curvature can be deduced formally from 
the exact expression for the disorder-averaged diffusion coefficient D (see below), which is 
valid in all dimensions d. The argument for the random-barrier model is more complicated. 
At the lowest temperatures, where the ratio of the width of the barrier energy distribution 
and the thermal energy is large, the diffusion coefficient is determined by the highest barrier 
along a critical path of bond percolation 0. The path may be constructed by selecting 
successively bonds with the lowest barriers possible until the "infinite" cluster in the lattice 
appears. If the temperature is raised, additional paths contribute to the diffusivity. These 
paths comprise also higher barriers, hence the apparent activation energy is increased. 

In the following section we describe the derivation of the exact result for D in disordered 
linear chains from an expression for the mean first-passage time. In Sec. 3 we give an 
analytical treatment of the combination of independent random barriers and random traps. 
Sec. 4 contains the EMA for the combined model in d > 3. In Sec. 5 we present our 
conclusions from the results. 

II. EXACT EXPRESSION FOR DIFFUSION COEFFICIENT IN D = 1 

Recently an exact expression for the asymptotic diffusion coefficient of single particles 
on linear chains with disordered transition rates has been derived pj-Q. In [0 a first- 
passage time method was used while in the mobility of particles on chains with periodic 
boundary conditions was derived from the linear response to a driving force. We follow 
reference [0] for the derivation. The basis is an exact expression for the mean first-passage 
time of a particle from site to site iV on a segment of a disordered chain Site is 
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regarded as reflecting, while site is an absorbing site. The expression reads 

N-l 1 N-2 1 N-l i -p 

A,.=0 ^ k,k+l jt=o ^ '=^.'=+1 i=k+l j=k+l 

Here Fj^j+i is the transition rate from site i to site i + 1; the transitions being restricted 
to nearest neighbors. Note that (P gives ioN for a particular (quenched) realization of the 
disordered segment. 

We now invoke the condition of detailed balance between two neighbor sites, 



PiTij = pjVji. (2) 

The thermal occupation factors pi are defined by 

_ exp{-(3Ei) 
{exp{-f3E,)} 

where Ei is the energy of site i, counted negative from a common origin. The curly brackets 
in d^) designate the average over the disordered site energies, for finite segments 

1 ^ 

{exp(-/3E,)} = ^exp(-/5E0. (4) 

1=1 

The occupation factors are proportional to the thermal equilibrium occupation probabilities 
of the sites; due to the normalization used they can be larger or smaller than one. 

The condition of detailed balance holds in equilibrium, and the occupation factors pi exist 
for finite segments when all nearest-neighbor transition rates Tij ^ 0. We have to require 
that unique occupation factors also exist in the limit ^ oo. This requirement excludes 
certain interesting models, for instance the Sinai model 0, from the further derivations. 
However, it is fulfilled for the models considered in the next section. 

When the condition of detailed balance (0) is introduced into expression (0) a consider- 
able simplification is achieved: 

N-l ^ N-2 N-l 

ioN = E r + E Pfc E • (5) 

k=0 fc=0 i=k+l 

Taking the disorder average of the equation, we find as the leading term in the limit of long 
segments, >> 1, 



This equation can be interpreted as the inverted relation between time and mean-squared 
displacements of random walks on disordered lattices, t = {2D)^^{X'^}. We hence deduce 
the following asymptotic diffusion coefficient from 

Under the assumption of the existence of unique occupation factors one can also derive, for 
large A^, using again detailed balance 

{Pon} - {ioNV oc N\ (8) 

Since {toiy}^ ~ A^'^, the relative dispersion of the disorder average of the mean first-passage 
times vanishes as N~^^'^ for long segments. One can say that {Iqn} becomes "sharp" in the 
asymptotic limit of large A^. In this way the use of the inverted relation (|^) to deduce the 
diffusion coefficient can be justified. 

The physical significance of the result (|^) is that the diffusion coefficient follows from 
thermally weighted transition rates. A numerical verification of (|^ for the Miller- Abrahams 



model PI in (i = 1 has been given in reference |10 



III. COMBINATION OF RANDOM BARRIERS AND RANDOM TRAPS: 

ANALYTICAL RESULTS 

A. Model and exact result in d = 1 

We now investigate a model for hopping diffusion of particles in lattices with combina- 
tions of random barriers and random traps. The random-barrier (RB) model is defined with 
symmetric transition rates Tij between neighbor sites, Tij = Tji. An Arrhenius form for the 
transition rates is assumed, 

rij = roexp{-f3Eij). (9) 



The energetic barrier Eij between sites i and j is a random variable and it is taken from a 
common distribution vb{E). To avoid problems with negative barriers we restrict the range 
of the barrier energies to > 0. The random-trap (RT) model is defined by rates Fj that 
originate from the sites i and are independent from the final sites. Also here the Arrhenius 
form is assumed, 

ri = roexp(/3E,) (10) 

where Ei is the energy of site i. The site energies will be counted negative and the range 
be restricted to E < 0. Again the individual site energies are selected from a common 
distribution ut{E). 

We introduce a combination of the RT and RB models by specifying the transition rates 
between two neighbor sites as 

r,, = roexp[-/3(E,,— Ei)] (11) 

The energy Eij with two site indices refers to the barrier i j while the energy E^ with one 
index refers to the site i. (To distinguish between barrier and site energies, we sometimes 
keep dummy site indices). A pictorial representation of the model is given in figure 1. Note 
that the model is different from the Miller- Abrahams model where the site energy of the 
terminal site j appears explicitly. 

The occupation factors pi are required in the weighted transition rates that appear in 
(0); they are given in (^. The weighted transition rates are then 

_ Toexp{-(3E,,) 

The numerator contains only barrier energies while the site energies only appear in the 
denominator, which can be evaluated directly. An immediate consequence of (|1^) is the 
expression for the diffusion coefficient of the RT model, 

Drt = ro{exp(-/?E,)}-i (13) 
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when barrier disorder is absent. Equation (|T3]) follows from (|^) in d = 1; it is valid in all 



dimensions [jTl|. Another consequence of Eq. (|T3|) is that the derivative of ln(D/jT) with 
respect to the inverse temperature (3 is given by the mean thermal energy. Since the mean 
thermal energy can only decrease with decreasing temperature, the slope of \n.[D ftril^)] is 
decreasing with increasing j3. 

An exact result is obtained from for the diffusion coefficient of the combined RB and 
RT model in = 1, 

Dcomb = To{exp{-pE,)}-'{exp{PE,,)}-\ (14) 
Equation (|1^ can be cast into another form, 

Dcomb = TtDrtDrb- (15) 
i 

The diffusion coefficient Drt exhibits generally downward curvature in an Arrhenius plot. 
Since the average that determines the diffusion coefficient Djib in ci = 1 has the same form 
as the one determining Drt-, downward curvature is also present in the Arrhenius plot of 
this coefficient. Also the diffusion coefficient of the combined model has then downward 
curvature m. d = 1. In other words, a compensation of the effects of random barriers and 
random traps is never possible in ci = 1. 



B. Results in higher dimensions 

The diffusion coefficient of particles in the random-barrier model in simple-square lattices 
is exactly known for symmetric barrier-energy distributions vsiE) [l^. For symmetric 
energy distributions the diffusion coefficient is 

DRB = T^eM-l3E) (16) 

where E is the medium value of the energy distribution. Hence Drb does not show curvature 
in an Arrhenius plot. Consequently, the diffusion coefficient of the combined model will 
exhibit downward curvature in simple-square lattices. 
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Approximations are necessary to derive the asymptotic diffusion coefficients in higher- 
dimensional disordered lattices. It is very plausible that the weighted transition rates (|I2|) 
should also be used in approximate treatments. Site-energy and barrier disorder can be 
treated independently when they are uncorrelated. Under this assumption, the numerator 
of Eq. (JT^) only contains uncorrelated barrier energies. The EMA gives reasonably accurate 
results for the RB model when the disorder is not very strong. Hence we will use this 
approximation to deal with the barrier disorder ind > 1. In the EMA, an effective transition 
rate Fe// is determined from a self-consistency condition [0,0, 

The transition rate F is a random variable, here it is taken according to Eq. ([1^). The 
rate F is symmetric as a consequence of detailed balance; this symmetry is required for the 
application of Eq. ([T7|). The curly brackets indicate the disorder average which extends over 
the barrier disorder. The diffusion coefficient is identical to Fe// and the lattice constant in 
the hypercubic lattices that we study is set unity. 

Since the exact result ( ]TB| ) for Djit is contained as a factor in the weighted transition 
rates, we have in EMA 

-L 

The EMA becomes exact in the limit of coordination number z —>■ oo, i. e., in infinite- 
dimensional disordered lattices. We obtain in this limit 

Dcomb = ro{exp{-pE,)}-\exp{-(3Eij)}. (19) 

The second average is proportional to the average over the rates corresponding to the random 
barriers. For the RB model and in ^ oo we have Drb = {^rb}, hence the product form 
of Dcomb (|15[) is also valid in infinite dimensions. We conjecture that (|T5|) is generally valid 
if the site and barrier energies are uncorrelated. Very recently Limoge and Bocquet incor- 
porated the backward correlations of the RB model and derived a self-consistency condition 
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for the effective transition rate. However, their self-consistency condition is different from 
the standard form ( p!7D of the EMA and does neither reproduce the exact one-dimensional 
(section [111 A] ) nor the exact two-dimensional result 



C. Possibility of complete compensation 



In this subsection we study the possibility for a complete compensation of the curva- 
tures resulting from the RT and RB models. For this purpose we consider first the case of 
infinite dimensions, d ^ oo. To obtain complete compensation of the curvatures, a special 
relationship between the site and barrier energy distributions must exist. To determine this 
relation, we require that the following equation represents simple Arrhenian behavior of the 
diffusion coefficient. 



with a temperature- independent activation energy Ecomb > 0. 

Let us assume that vt{E) is restricted to the energy interval [— -Ec, 0] with a positive 
quantity. If we identify E^. with Ecomb we can show by simple manipulations of the left-hand 
side of Eq. (^) that it is satisfied when 



The barrier distribution is then restricted to the interval [0, Ecomb] and it is simply the trap 
distribution shifted by Ecomb- The argument implies that complete compensation, in d ^ oo, 
is only possible for energy distributions that are restricted to a finite interval. Of particular 
interest are distributions that are symmetric about the midpoint of the interval. Of course 
it is not necessary that the distributions are unequal zero on the whole interval. 

For finite dimensions 1 < d < oo the diffusion coefficient Dub has to be determined 
by approximations, for instance by the EMA. It is a complicated functional of the energy 
distribution ub{E) and no general result can be obtained. It seems very implausible that 



ro{exp(-/3£'i)} ^{exp(-^£'jj)} = Fq exp(-/3£'com6) 



(20) 




(21) 
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complete compensation can be achieved in general. However, it appears always possible, in 
dimensions d > 3, to choose the energy distributions vxiE) and vsiE) in such a way that 
an approximate Arrhenian behavior is achieved in a restricted temperature interval. 



IV. COMPARISON OF EMA RESULTS AND NUMERICAL SIMULATIONS 

In this section we present EMA results for the combined model in c? = 3 and 5 for different 
forms of the disorder and various temperatures. The results are compared with Monte-Carlo 
simulations of the diffusion coefficient by measuring the mean-square displacement 
where the asymptotic diffusion coefficient D is given by 

{r^it)] ^2dDt. (22) 

Actually we present the whole proportionality factor 2dD as the diffusion coefficient in the 
figures to make results in different dimensions comparable. Compensatory effects between 
random barriers and random traps are now discussed. 

We first calculate DfJ^^ for uniform distributions of energies ( for the RB model see 

H), 

d—l 2 smh(2^pcrB) exp^paT) — 1 



where 



^ -gt <E<0 (RT) 
u{E) = \ ^ 0<E<aB (RB) 



otherwise 

By comparison with the numerical simulations we find that using 

d — l , 
ctt ~ ~~^^B (24) 

approximate Arrhenian behavior is reached. In figure 2 results of Monte Carlo simulations 
are shown together with the EMA results (p3|) with parameters chosen according to 
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The upward curvature for the RB model is stronger in d = 5 than in c? = 3 although the width 
of the distribution is the same. The compensation works in the case of a stronger curvature 
for the RB and RT models because for larger coordination numbers compensation becomes 
generally easier. The effect of the coordination number on the possibility of compensation 
of barrier and site-energy disorder was already noted in . 

We consider as a second energy distribution the Gaussian distribution, 



u(E) = exp (-— — ^ 1 with E { 



(25) 



< 0, a = aT for RT 
> 0, a = (Tb for RB 

We have to cut off the tails of one side of the distributions to get only negative respectively 
positive energy values. We decided to cut off 5% of the distributions. Therefore we take E 
according to 

E = -(Terf-^(0.95) (26) 



where erf() is the error function. For this distributions the self-consistency condition (|T] 
has to be solved numerically. Results are shown in figure 3. For this example the two curves 
for the RT model lie on top of each other because we took the same width in ci = 3 and 
5 and then determined the width of the barrier distribution to get compensation in the 
combined model. Again the upward curvature of the RB model is stronger in d = 5 than in 
d = 3 although we took a smaller width of the barrier distribution in the higher dimension. 
Nonetheless compensation of the curvature is still possible and the numerical procedure to 
find the correct parameters of the distribution for the best possible compensation is even 
more stable in the higher dimension. Though the curves of the combined model in figure 3 
look linear, the numerical results show that no perfect compensation is possible. 



V. CONCLUSION 



Our theory of the diffusion coefficients of models with combined site-energy and barrier 
disorder was based on weighted transition rates where the thermal site occupancies enter 
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explicitly. Due to the exponential form of the assumed Arrhenius law for the individual 
transition rates, the contributions from the site energies and the barrier energies factorize. 
If the site-energy and barrier disorder are uncorrelated, the resulting diffusion coefficient 
factorizes into the random-trap and random-barrier contributions. A consequence is that 
compensation of the effects of random site-energies and random barriers on the curvature in 
an Arrhenius plot of Dcomb versus f3 is not possible in d = 1 and d = 2. The simple-square 
lattice represents a boundary case. 

Partial compensation is possible in a finite temperature interval in 3 and higher dimen- 
sions, if the strength of the disordered site and barrier energies is adjusted properly (stronger 
disorder in the barriers than in the site-energies for finite dimensions d > 3). 

In this context we should mention that our conclusions refer to the hypercubic lattices 
(square, simple-cubic, etc.). The relevant parameter for the temperature dependence of the 
diffusion coefficient in the RB model is the coordination number z. Lattices with z ^ 2d 
were investigated by Limoge and Bocquet 0. 

We conclude that partial compensation of the effects of random barriers and random site 
energies is possible in c? = 3, in a finite temperature interval, if two premises are fulfilled, 

i) Assumption of independent site energies and barrier heights, 

ii) Properly adjusted strength of the disorder. 

Point i) may be somewhat alleviated by including short-range correlations, but point ii) 
seems to be generally necessary. We leave the question open whether the points i) and ii) 
are reasonable descriptions of real amorphous substances. 

Instead we point to one serious deficiency of the present theories of diffusion in disordered 
system. The present theories are based on regular lattices with the disorder put into the 
transition rates. Real materials have topologically disordered structures. There exist a range 
of coordination numbers of the equilibrium sites, and the jump distances vary considerably. 
The treatment of these effects remains as a task for the future. 
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FIGURES 



FIG. 1. Schematic representation of the combination of random barriers and random traps. 

FIG. 2. Diffusion coefficient in the RT, RB, and combined model, for uniform distributions of 
the energies. The parameter Fq is set unity. The different symbols represent simulation results for 
the RB (+), RT (o), and the combined model (□) with ut = 3.0,3.2 for d = 3 (full curves) and 
d = 5 (dashed curves) respectively and ctb = 4.0 in both cases. The curves represent the EMA 
result for 2dDcomb (HI)- 

FIG. 3. Diffusion coefficient in the RT, RB, and combined model, for the Gaussian distributions 
of energies. The parameter Tq is set unity. The different symbols represent simulation results for 
the RB (+), RT (o), and the combined model (□) with as = 4.0, E ~ 6.58 for d = 3 (full curves), 
as ~ 1.74, E ~ 2.86 for d = 5 (dashed curves), and ax = 1.2, E ~ —1.97 in both cases. The 
curves represent the EMA result for 2dDcomb (p3|). 
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